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Abstract 

We present a model of the Galactic Habitable Zone (GHZ), described in terms of 
the spatial and temporal dimensions of the Galaxy that may favour the development 
of complex life. The Milky Way galaxy is modelled using a computational approach 
by populating stars and their planetary systems on an individual basis using Monte- 
Carlo methods. We begin with well-established properties of the disk of the Milky 
Way, such as the stellar number density distribution, the initial mass function, the 
star formation history, and the metallicity gradient as a function of radial position and 
time. We vary some of these properties, creating four models to test the sensitivity 
of our assumptions. To assess habitability on the Galactic scale, we model supernova 
rates, planet formation, and the time required for complex life to evolve. Our study 
improves on other literature on the GHZ by populating stars on an individual basis and 
by modelling SNII and SNIa sterilizations by selecting their progenitors from within 
this preexisting stellar population. Furthermore, we consider habitability on tidally 
locked and non-tidally locked planets separately, and study habitability as a function 
of height above and below the Galactic midplane. In the model that most accurately 
reproduces the properties of the Galaxy, the results indicate that an individual SNIa 
is ~5.6x more lethal than an individual SNII on average. In addition, we predict that 
~1.2% of all stars host a planet that may have been capable of supporting complex 
life at some point in the history of the Galaxy. Of those stars with a habitable planet, 
~75% of planets are predicted to be in a tidally locked configuration with their host 
star. The majority of these planets that may support complex life are found towards 
the inner Galaxy, distributed within, and significantly above and below, the Galactic 
midplane. 

Keywords: astrobiology, Galactic evolution, Galactic Habitable Zone, habitability, su- 
pernovae 



1 Introduction 



Studies of habitability on the Galactic scale are gaining attention as extrasolar planetary 
searches and advanced models of the evolution of the Milky Way galaxy improve our under- 
standing of the prerequisites for life. The GHZ is the area in the Galaxy that may favour 
the emergence of complex life. The inner boundary of the GHZ is thought to be defined by 
hazards to planetary biospheres, and the o uter boundary is set b y the minimum amount of 
metallicity required for planet formation ( Gonzalez et al.l 120011 ). Given these factors, the 
GHZ is constrained by the morphology, stellar populations, and chemical evolution of the 
Milky Way. The GHZ is defined as a concept analogous to the Circumstellar Habitable 
Zone (HZ), which is the region around a star where water can remain in a liquid state on a 
planet's surface. 

To assess habitability on the Galactic scale, we model planet formation, supernova rates, 
and the time required for complex life to evolve. Estimates in predicting the number of hab- 
itable planets in the Galaxy are in their infanc y at present . Can didate habit able p lanets 
have been identified in the Gliese 581 syst em (Selsis et ah . 20071 : iVogt ct aU |2010| ). and 
in the Kepler list of candidate exoplanets ( Borucki et al.l 120111) . Models can provide in- 
sight into the number and distribution of habitable planets that are likely to exist in the 
Galaxy. Prior studies suggest that there is a cor r elation betwe e n high-metallicity environ- 
ments and planet formation J Ida and Lin j_ 2004 ■ ISantos et al.l 2004 : iFischer and Valenti , 
I2OO5I: iGrether fc Lineweaveii 120071: ISozzetti et al 1 . 120091 : IJohnson et all l2010t ). The meteJ- 
licity gradient found in the disk of the Galaxy has been measured, and suggests that the 



number of planets depends on radial position. Regions greatly affected by the detrimental 
impacts of supernovae (SNe) are defined by proximity to these transient events. There- 
fore, the stellar density and star formation history (SFH) directly influence the frequency at 
which a planet's biosphere is negatively affected. Finally, in combination with the periods 
that SNe are expected to sterilize life in surrounding solar systems, enough time must elapse 
such that complex life can e volve. 

Research on the GHZ bv iLineweaver et al. (120041) sug gests that the GHZ is an annular 
region that expands with time. However. IPrantzosI ( 2008 ) argues that the entire disk of the 
Galaxy may be suitable for complex life. From our work we predict that at the present 
time, the entire disk of the Galaxy is hospitable to complex life; however, the regions that 
permit the majority of complex life are contrary to other predictions in the field. 

Predicting where complex life might be favoured in the Galaxy follows from making 
assumptions about the prerequisites for life. Whether or not life exists beyond the context 
of our solar system cannot be validated at present. Our definitions of habitability in the 
Galaxy do not imply any claims of inhabitance. 

We model the GHZ with the factors described above. We begin in § [U by outlining 
the construction of our model Milky Way galaxy. In § [3] we describe the factors that 
influence habitability on the Galactic scale. These factors include: the dangers caused by 
SN sterilizations, how the metallicity gradient influences planet formation, and the timescale 
requirements for the emergence of biologically complex life. Results and a comparison with 
other studies are discussed in §|31 Finally, in §[5] we present the conclusions drawn from our 
study of habitability in the Galaxy. 



2 A Model of the Milky Way Galaxy 

In this section we describe our model of the Milky Way. We assign the birth date and 
metallicity of each star in a se l f-cons istent way by matching the SFH and the metallicity 
gradient of iNaab and Ostrikeij (|200g) . To assess the sensitivity of our model, we vary the 
initial mass function (IMF) and the stellar number density distributions. Varying the stellar 
number density distributions creates differing profiles of the star f ormation rate (SFR) within 
the Galactic disk. We do not try to match the number density of lNaab and Ostriken (|2006 ) , 
as their model does not describ e the vertical component of th e disk. Rather, we match the 
number density distributions of iGarroll and Ostlid ( 20061) and lJuric et al.l ( 2008 ). to permit 
us to study the Galaxy as a function of distance above the midplane and radial distance 
from the Galactic centre. This approach yields a three-dimensional model that allows us to 
study habitability as a function of radial distance and height above the midplane. 



2.1 Stellar Mass and Main Sequence Lifetime 

This subsection outlines how a mass and main sequence lifetime are assigned to each star 
in our model Milky Way. We assign stellar masses using a Monte Carlo technique such that 
the r esulting distribution follows a power-law and matches an IMF with a = 2.35 (jSalpeter . 
19551 ). The maximum and minimu m stellar masses for main sequence stars are defined as 



100 Mq and 0.08 Mg respectively (JKroupa et al.l . ll993l) 



This IMF is consistent w i th ob servations of star-forming regions; h owever, there are 
uncert a inties (Leitherer et al. . 19991). and oth er IMFs have been proposed. iMiller and Scald 
( 19791 ) ■ iKroupal (2001), and IChabrieiJ ( 2003 ) suggest that at subsolar masses, a. should be 



lower than 2.35, flattening the IMF. iNaab and Ostrikeij ( 2006[) test both the Salpeter and 
Chabrier IMF (flattening at subsolar masses) in their model of the Milky Way and find that 
both are in agreement with observations. One of the IMFs we uti l ize is the Salpeter IMF 
to remain consistent with the model Galaxy bv lNaab and Ostrikeij ( 20061 ). 

To test the impact of the IMF flatteni ng, in comparis on to the Salpeter IMF at subsolar 
masses, we implement the IMF given bv iKroupal ( 2001 ). which is described as a two part 
power law function. The value a = 1.3 when 0.08 < M < 0.5, and a = 2.3 when M > 
0.5. The steep slope at subsolar masses in the Salpeter IMF suggests that the Galaxy is 
populated with many more low mass stars than the Kroupa IMF. ^^^^^^_^^^^^^^ 

The mass corresponds to a main sequence lifetime defined by [Hansen and Kawaleij 



(119941) 
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where ?7Zf7) = 1 is the Sun's mass, T^q = 11 is the Sun's main sequence lifetime in Gyr 



(jSackmann et al.l . ll993i ). and m is the star's mass determined from the Salpeter or Kroupa 



IMF. 



2.2 Star Formation History 

This su bsection outlines when star s form in our model of the Milky Way. We adopt the SFH 
given bv lNaab and Ostrikeij ( 2006[) . who model the SFH in the Galactic disk over time. Their 
model yields an inside-out formation of the Galaxy, with an early burst of star formation in 
the in ner Galaxy that declines over time. Close inspection of Figure 9 in lNaab and Ostrikei 
(|2006[ ) reveals that the re are fluctuations in t he SFH, as natural processes such as the 
passage of spiral waves ( Hernandez et all . l2000f ) and close encounters with the Magellanic 
Clouds (JRocha-Pinto et al.l . 120001 ) trigger star formation. These sophisticated mechanisms 
that cause fluctuations in star formation are not modelled in the presen t study. Fitting the 
curves and interpolating between them in Figure 6 (SFR vs. radius) from lNaab and Ostrikei 
(|2006l ) establis hes when the stars form in our model. W e recognize that other SFHs have 
been reported (JAumer and Binnevj . I2009J : IFu et al.l . 120091) , although we do not model them 
in this study. The SFH adopted here varies by radial position (R), but is independent of 
the vertical (z) position in the Galaxy. 



2.3 Distribution of Stars Within the Model Milky Way Galaxy 

We utilize the distribution of stars found bv lCarroll and Ostliel ()2006l ) and ljuric et al.l (|2008h 
to determine where the stars form in our model of the Milky Way. We vary the number 
density distributions to test how sensitive the model is to differing profiles of this parameter. 



2.3.1 Distribution of Stars bv ICarroU and OstlQ (|2006l ) 



Their formula describes the number density of stars as a function of distance above the 
midplane and radial distance from the Galactic centre, as follows: 



n(^z,R) = no(e"^/^""" + 0.085e~^/^"'"'=)e"-"/''«. 



(2) 



The quantity n is the number of stars per unit volume (pc^). The coordinate z is the vertical 
height above the midplane of the Galaxy, and R is the radial distance from the Galactic 
centre. The constant Kr = 2.25 kpc denotes the radial disk scale length. The thin disk lies 



within the central plane of the Galactic disk with a scale height of zthin — 350 pc, and the 
thick disk has a vertical scale height of zttick —1000 pc. The quantity no (number of stars 
per pc^) is the normalization that we will select below. 

We normalize the distribution of st ars based on an estimate of th e total disk mass in the 
Milky Way. The disk mass estimate of iBinnev and Tremaind ( 20081 ) yields 4.2 x IO^^Mq in 
disk stars (0.3 x lO^^M© has been subtracted, as it corresponds to gas mass). We are able 
to match this disk mass using our Salpeter IMF described in ij 12. II with a value of rig =5. 502 
stars pc~'^. This normalization yields ~150 billion stars in the disk of the Galaxy. With 
regards to the Kroupa IMF, we normalize no=1.957 stars pc~^, yielding ^^50 billion stars 
in the disk of the Galaxy. 

To ensure our distribution is appropriate, we compare our estimates to the number 
density of the local neighbourhood. For the Salpeter IMF, t his normalization overpre- 
dicts by a factor of ^1.4 the local number density reported bv iReid et al.l (|2002l) . Other 



Carroll and Ostlid . 120061 ) are even higher than that of 



disk mass estimates ( f or exa mple, 

Binnev and Tremaind ([2003), and would overpredict the local density to a greater degree 



than the normalization utilized here. Note that the 40% discrepancy refers only to the local 
neighbourhood, and would be unacceptable if it applied t o the entire disk. For the Kroupa 
IMF, the local number density is ~50% of that found by iReid et al.l (120021). The Salpeter 
and K roupa IMFs, paired with the stellar number density distribution of ICarroU and Ostlie 
(|2006l ) are referred to as Model 1 and Model 2 respectively. 



2.3.2 Distribution of Stars bv lJuric et al.l (|2008h 



Their formula describes the number density of stars as a function of distance above the 
midplane and radial distance from the Galactic centre, as follows: 



Pd{R. Z) = pd{R, Z; Li, H^) + fpoiR, Z; L2, H2 



(3) 



where 



PDiR,Z;L,H) = pDiRe,0)e 



Rq/L 



X e^ 



-^). (4) 

The quantity pn is the number of stars per unit volume (pc'^). The coordinate Z is the 
vertical height above the midplane of the Galaxy, and R is the radial distance from the 
Galactic centre. We utilize the values HI = 300 pc, LI = 2600 pc, H2 = 900 pc, L2 = 
3600 pc, and / = 0.12, corresponding to the thin disk scale height and length, the thick 
disk scale height and length and the thick-to-thin disk density normalization. 

Using the Salpeter IMF, we are able to match our disk mass (4.2 x lO^^M©) by normaliz- 
ing pd{Ro,0)— 0.237 stars pc""^. We normahze p_d(_R0, 0)=0.084 stars pc"'^, corresponding 
to the Kroupa IMF. 

The norma l izatio n for the Salpeter IMF overpredicts the local number density reported 
bv iReid et al.l ( 20021 ) by a factor of ^2. When implementin g the Kroupa IMF, the local 
number density is ~70% of that found bv iReid et al.l (|2002). The Sa l peter and Kroupa 
IMFs, paired with the stellar number density distribution of IJuric et al.l (|2008t ) are referred 
to as Model 3 and Model 4 respectively. 



2.3.3 The R and z Dependence of the Number Density Distributions 

Our model Galaxy has both R and z dimensions. The vertical motions of stars passing 
through the midplane, which are not included in this study, would serve to reduce (but 



not remove) trends in the z dimension; however, we do not expect trends in habitabihty 
to be significantly diminished. Stars spend different amounts of time above and below the 
midplane, as all stars do not bob up and down to the same extent. For example, the Su n 
is expected to remain within 70 pc above or below the midplane ( Gies fc Helsell . |2005[) . 
Therefore, the Sun is well within a scale height of the midplane throughout its vertical 
motions. We do not believe that the absence of vertical and similarly, radial motions in our 
model to substantially affect trends as a function of z or R. 

Our model of the Milky Way ignores binary star systems, spiral arms, the orbits of 
stars, the production of stars in clusters, and the accretion of satellite galaxies over time. 
These are observable properties of the Galaxy that could have an impact on habitabihty, 
but are not implemented in the present study. While other studies use more sophisticated 
simulations of the Milky Way, we populate stars on an individual basis. Our approach is 
reasonable, since it reproduces the global properties of the Galaxy. 

The SFR produced in our model as a function of radial distance and time is shown in 
the upper (Model 1) and middle (Model 4) panels of Figure [1] As a result of the inside-out 
nature of Galaxy formation, the inner Galaxy has an early burst of star formation that 
decreases with time. The SFR around the solar neighbourhood has remained constant over 
the past few Gyr, and the SFR has been increasing in the outskirts. 



2.4 Metallicity 

Through the continual births and deaths of generations of stars, the abundance of metals 
increases in the Galaxy over time. Metals are those elements heavier than helium, and 
are important as they are the building blocks of terrestrial planets. Stars that form at 
a later date have a greater chance of having terrestrial planets than earlier generations 
of stars. We use metallicity as a proxy for planet formation , which is consistent with 
other literature on the subject (Ida and Lin . 2004 : Santos et al.l. 2004t iFischer and Valentil . 



20051 : iGrether fc Lineweaveill2007l : ISozzetti et al.l . T2009Hjohnson et alJ . r2010[ ). A metallicity 
profile is used to assign an appropriate metallicity in terms of [Fe/H] abundance to each 
star in the model^. There are a number of models that map the evolutio n of the metallicity 
gradient in the Galaxy. We utilize the metallicity profile in Figure 11 of iNaab and Ostriker 
(12006,). We use a fixed metallicity for stars born at a given radial position and time, which 
is consistent with lNaab and Ostrikeij ( 20061 ). The average metallicity of the stars produced 
in terms of radial position and time in our model is shown in the lower panel of Figure [TJ 
and the distribution of metallicities for all stars that form at the solar radius (R~8 kpc) 
is shown in Figure [5] The model bv lNaab and Ostrikeij (2006) is consistent with literature 
regardin g the metallicity di stribution of stars within the solar neighbourhood. At this 
location. iRobles et al.l ( 20081 ) find that 35% of local stars have Z>Zo, whereas 41% of the 
stars in our model have Z>Z0. Metallicity varies as a function of K and time in our model, 
whereas a more complex model could also depend on z. Metallicity trends as a function of 
z may be minimal, as studies suggest that there is little variation in the vertical gradient 
within the thick disk ( Gilmore et al.l . Il995l) . 

The chemical evolution of the Galaxy is affected by radial mixing. Radial mixing has 
the effect of modestly flattening the metallicity gradient across the Galactic disk; however, 
the gradient is steep enough such t hat planet formation is not iceably varied throughout the 
disk of the Galaxy. Therefore, like iNaab and Ostrikeij ( 20061) we do not account for radial 



^If a star is formed before 2 Gyr, it is assigned the metallicity of a star that formed at 2 Gyr. 



mixing in our model, as this mechanism would not be responsible for significantly changing 
the observed trends in the results. 



2.5 Varying Models of the Milky Way 

We vary the number density distribution and IMF to see how sensitive our model is to 
different estimates of the properties of the Milky Way. Table [T] gives an overview of our 
models of the Galaxy. We focus on Models 1 and 4 in our results, as they ar e the most 
consis tent with the number density of the solar neighbourhood as reported by iReid et al. 
(J2002i ). Model 1 overpredicts the local stellar number density by a factor of ~1.4, and 
M odel 4 unde r estima tes the local stellar number density, and is ~70% of the value found 



bv iReid et al.l ( 2002 ). These two models bracket the stellar number density in the solar 
neighbourhood. 



2.6 Model Implementation and Properties of the Galaxy 

Our model of the Milky Way disk has a radius of 15 kpc. We present our findings for the 2.5- 
15 kpc region of the disk, as the bulge overlaps the disk at R <2.5 kpc. We take advantage of 
the azimuthal symmetry within the Milky Way by modelling only a 1° sector of the disk. To 
eliminate uncertainties with this approach we introduce periodic boundary conditions when 
processing SNe sterilizations. Stars are processed on an individual basis; when processing 
SNe sterilizations, we keep track of the times in which SNe are sufficiently close to each star. 
We iterate through all of the stars in the model without using a brute-force approach. 

3 A Model of Habitability Within the Milky Way Galaxy 

Various properties were reviewed in §[5] that are used to construct a model of the Milky Way 
for the purposes of assessing habitability on the Galactic scale. In the following subsections 
we outline some prerequisites and conditions for land-based animal life. These prerequisites 
are used to assess if there is a position and time in the Galaxy that favours the emergence 
of complex life. 



3.1 The Effects of Supernovae on a Planet's Biosphere 

Cosmic rays, 7-rays and X-rays emitted by a supernova can lead to the sterilization of 
life in surrounding solar systems (for an overview of astrophysical radiation sources see 
Melott fc Thomasl . I2OIII) . The depletion of ozone in a planetary atmosphere caused by a 
SN exposes a planet's surface to more radiation from its host star ( Gehrels et al.U 2003l). We 
model Type II SNe (SNII) and Type la SNe (SNIa) independently to reflect the differences 
in their formation rates, luminosities and consequent sterilization distances. SNe that are 
not SNIa or SNII are ignored, as they are respons i ble fo r only a small fraction of SNe 



that oc cur in all types of galaxies (ICappellaro et al. , 119991 ) . This method differs from the 
work of iLineweaver et al.l ( 20041 ) and IPrantzosI ( 20081 ). Our approach treats both SN types 
separately and permits a range of sterilization distances for SNII and SNIa. We model 
SNe in a self-consiste nt way by selecting SN pr ogeni tors from within our pre-existing stellar 
distribution, whereas ILineweaver et al.l (200J) and IPrantzod ( 20081) use a time-integrated 
supernova danger factor which depends on the supernova rate, but is normalized in an 
arbitrary fashion. 



3.1.1 Supernovae Sterilizations in Other Studies 

The previous work of iLineweaver et al.l ( 2004 ) and IPrantzoa ( 2008 ) use the SNe rate in the 



Galaxy to determine the probabihty that hfe on a planet survives SNe explosions across 
the Galactic disk. They normalize the danger posed by SNe to the time integrated SN rate 
of the solar neighbourhood. They suggest that if the time integrated rate is 4x that of 
the solar neighbourhood, the probability of surviving SNe explosions is 0, and if the time 
integrated rate is <0.5 that of the solar neighbourhood, the probability of surviving SNe 
explosions is 1. 

Whether a planet survives a SNe explosion is linked to the SNe rate; however, the time 
in which a SN occurs in a planet's history is also a major factor. For example, if a planet 
is irradiated by multiple SN explosions early in its history, the sterilizations should not 
preclude it from being habitable at later dates. In an analogous scenario, the late heavy 
bombardment has not prevented the Earth from attaining complex life. We test to see when 
each individual planet was irradiated by SN events, which permits us to more accurately 
determine if a planet has the potential to host complex life. This is described in Section [2111 
Furthermore, we do not model a SN rate; however, if we were to test the rate at various 
positions, ours would depend on R and z, whereas the related work is only dependant on 
R. 



3.1.2 Supernovae Sterilizations in the Present Study 

Stars wit h a main s e quenc e mass greater than wS Mq become SNII at the end of their 



lifetimes ()KennicuttJ . 119841) . In our model, all stars with a main sequence mass >8 Mq 
become SNII after their respective main sequence lifetime expires. We do not account for 
the delay between the end of the main sequence lifetime of a star and the time before it 
explodes as a SN; however, the time delay before a SN occurs is relatively small, as a large 
majority of a star's lifespan is spent on the main sequence. Models 1 and 3 (Salpeter IMF) 
yield a Galactic SNII rate over the past Gyr of 1.98-2. 40x 10^^ yr^^, consistent with the SNII 
rate of ICappellaro et al.l (|l993l ). who find a Galactic SNII rate of 0.4-2xlO"Vi""^- Models 
2 and 4 (Kroupa IMF) yield a Galactic SNII rate ov er the past Gyr of 3 70-4 40x10'^ yr^^. 



consistent w ith the SNII rate by a factor of ~2 of ICappellaro et al.l ([19931). The average 



SNII rate in ICappellaro et alj (Il993[) is calcu lated bv using the luminosity of the Galaxy. 

Theoretical estimates bv iTutukov et al.l (|l992[ ) predict a SNII rate of 1.96-3.35x10^2 
yr^^, which is closer to the predicted values in the present study. The Kroupa IMF leads 
to more high mass stars than the Salpeter IMF; therefore, the Kroupa IMF yields more 
SNII. ICappellaro et al. (1993il note that six Milky Way SNe have been observed within the 
last millennium, but many additional Galactic SNe have gone undetected due to absorption 
by dust in the disk. Predictions of the true SNe rate vary considerably in the literature. 

The type of SNIa focused on in our model is the single degenerate case, wherein a white 
dwarf in a binary system grows in mass due to accretion from a binary companion. When 
the Chandrasekhar limit of l AMc^ is reached, the star explodes as a result of insufficient 
electron degeneracy pressure. iPritchet et al.l (|2008l ) estimate that ^1% of all white dwarfs 
become SNIa, independent of mass. In our model, all stars that are white dwarf candidates 
have main sequence masses of 0.08 Af© < M < 8 Mq, and are assigned a 1% chance of 
becoming SNI a after their main seq uence lifetime has expired. The detonation delay times 



implicit in the | Pritchet et al.l (j200 



model are consistent with the delay time estimates of 

3 ,,,.-1 



Raskin et ahl (|2009t ). Models 1 and 3 (Salpeter IMF) yield a SNIa ra te of 4.0x10'" yr~" ove r 
the past Gyr, consistent within a factor of ~ 1.5 of the SNIa rate of iMeng and Yana ()2010l ). 



who find a Galactic SNIa rate of 2.25-2.9x10"^ yr"^ Models 2 and 4 (Kroupa IMF) yield 
a SNIa rate of 6.6xl0~^ yr~ ^ over the past Gyr, consistent within a factor of ^2.5 of the 
SNIa rate of iMeng and YansJ ( 20101 ). The rate of SNIa is lower than that of SNII because 
the majority of SNIa candidates are very low mass stars and long timescales are required 
before they evolve off the mai n sequence. We r ecognize that the double degenerate channel 
is responsible for some SNIa (JMeng and Yana, 2010); however, obser vations indicate that 



the single degene r ate ch annel is dominant ([Parthasarathv et al.l . 120071 ) 



Gehrels et al.l ()2003() find that at a distance of < 8 pc, a SNII will deplete the ozone in 



a planet's atmosphere to the point where the UV flux received from the local host star will 
have a sterilizing effect on any land-in habiting life that exis ts on the planet. We assume that 
the sterilization distance outlined bv lGehrels et al.l ( 2003 ) refers to an average SNII and is 
just sufhcient to s terilize life within 8 pc. W e scale the sterilization distances of all SNII and 
SNIa accordingly. iRichardson et al.l (J2002I) address the distribution of absolute magnitudes 
for SNII. Considered in the study are 54 observed SN events in external galaxies, and the 
mean absolute magnitudes in the B band (Mb) are found for each of the SNII subclasses 
(SNII-L, SNII-P and SNIIn). We fit the data across the SNII subclasses to determine the 
distribution of absolute magnitud es of SNII. In add ition, we use the distribution of SNIa 
absolute magnitudes reported bv IWang et al.l (|2006l ). which is based on 109 SNIa with a 
mean Mb of -19.34. 

Given these distributions of SNe absolute magnitudes, we compute the sterilization dis- 
tances for each SNII and SNIa using the following: 



''SN 



8 pc X ^/lO-OA{MsN-M,tci)^ 



(5) 



where dsN is the resulting sterilization distance for a given SN, Msn is its absolute mag- 
nitude, and Mstd = —17.505, which is the absolute magnitude of the average SNII that we 
assume is just sufficient to sterilize life within 8 pc. 

Our resulting distribution of sterilization distances, corresponding to the distribution of 
absolute magnitudes of SNII and SNIa, are shown in Figure [S] SNIa are more luminous 
than SNII on average, and therefore are expected to sterilize planets at a greater distance 
than SNII. 



3.2 Metallicity and Planet Formation 

Predicting a region that may favour complex life in the Galaxy would be simplified given a 
census of habitable planets in the Milky Way. The Kepler mission will even tually yield an 
estimate of this nature. Of the 1202 Kepler exoplanet candidates discussed in lBorucki et al 



(J2011I ). 6 exist in the HZ that are less than twice the size of the Earth. In the meantime, we 
use the metallicity-planet cor r elatio n ( Ida and Linl.l2004t Santos et al.l. 2004 ;^ Fischer and Valenti 



2005t iGrether fc Lineweaveii l2007t ISozzetti et ail , l2009l : Ijohnson et all I2OI0I) in combina 
tion with our model of metallicity within the Galaxy and a model of solar system formation 
to predict where habitable planets will be found. This subsection describes how habitable 
planets are assigned to individual stars in our model Milky Way. 



3.2.1 Planet Formation in Other Studies 



Lineweaver et al 



(J2004l) calculate the probability of terrestrial p l anet formation in their 

model of the GHZ as influenced by iLineweaveii ( 20011 ). iLineweaven (2001j suggests that too 



little metallicity will inhibit the formation of planets and too much metallicity will produce 



giant planets that may hinder the formation or survival of habitable planets. Giant planets 
may reduce the number of habitable planets through the following processes: 1) the planet 
will be ejected out of the solar system or into its host star through the process of gravitational 
interaction, or 2) the protoplanetary disk will be consumed by the migrating giant planet, 
such that there is not enough ma terial to form ter restrial planets, or the migrating giant 
planet accretes habitable planets. iLineweaver et al.l ([2004) find that the potential of a star 
to harbour Earth-like planets increases linearly until a high metallicity abundance suggests 
the formation of massive planets. As the probability of having massive planets around a 
host star incre ases, the probab ility that Earth-like planets survive decreases. 

Converselv. iPrantzoa (J2008l) argues for a fixed planet formation probability of 40% above 
0.1 Zq. This is justified by the statement that Earth- like planets may be common in low- 
metallicity environments. The metall icity integrat ed probability of forming an Earth-like 



planet is the same as in the work by ILineweaver et al.l (|2004l ) . IPrantzosJ (|2008l ) estimates 
th e probability of form i ng Ho t Jupiters using the metallicity-planet correlation quantified 
bv iFischer and Valentil ( 2005 ). This work is also utilized in our model. The probability of 
having a star with a habitable planet without a Hot Jupiter that would de stroy it extends to 



lower and higher metallicity ranges than that of ILineweaver et al.l ( 20041 ). For comparative 



purposes, the probability of forming an Earth-mass planet in the model bv IPrantzod ( 2008[ ) 
is much higher than the probability adopted in our model. 



3.2.2 Planet Formation in the Present Study 

The probability of forming a habitable planet or a migrating massive planet (hereafter re- 
ferred to as a Hot Jupiter) is es timated using the m etallicity of each hos t star in the model. 
We use data from simulations bv llda and LinI ( 20051) as well as research bv lFischer and Valenti 
(J2005) to estimate the number of habitable planets and Hot Jupiters in our model. 

The probability of forming a gas giant planet around a. star is utilized to e s timat e the 
probability of forming a habitable planet or Hot Jupiter. iFischer and Valentil ( 20051 ) find 
that the probability of forming a gas giant planet is 



P{planet) = 0.03 x lO' 



2.0[Fe/H] 



(6) 



ISantos et al.l ( 20041 ) also quantify the planet- metallicity correlation, and find that a fiat 
metallicity tail exists, wherein the probability of forming a planet is constant (3%) at sub- 
solar metallicities ([Fe/H]<0). 



The Probability of Forming a Hot Jupiter 

We assume that the probability of forming a Hot Jupiter is 3% at sub-solar metalli cities and 
that t he probability of forming a Hot Jupiter at [Fe/H]>0 is given by Equation (6] IPrantzos 
(|2008f ) estimates the probability of forming a Hot Jupiter using a similar method. 



The Probability of Forming a Habitable Planet 

To link the metallicity-planet correlation of gas giant planets to habitable planets, we employ 
a model of solar system formation that contains both habitable planets and planets that 
w ould be d etectable by the radial velocity technique (short period gas giants). The model 
of llda and Lin (2005) produces short period gas giants and habi t able p lanets in roughly 
equal numbers. We note that in the work of IFischer and Valentil ( 20051 ). nearly all of the 



planets stu died satisfy the rad ial velocity criteria of llda and LirJ (|2004l ). In particular, in 
Figure 2 of llda and LinI (|2005r ) the distributions of the semimajor axis and mass of planets 
around stars of 0.2, 0.4, 0.6, 1.0 and 1.5 Mq are shown. Note that we assume a planet is 
considered habitable if it exists in the HZ around its host star (allowing liquid water on the 
planetary surface), a nd has 0.1 M^ < Mp < 10 M^, where Mp is the mass of the planet 
( Kasting et al.l Il993l ). A planet i s within the HZ if it s semimajor axis (a) lies within the 
range 0.8 - 1.5 x {M^/MqY AU (JKasting et all . Il993l ). where M* is the mass of the star. 



To associate the metallicity-planet correlation for gas giant planets to habitable planets, 
we find the ratio sof HZ planets wi th 0.1 Mg < Mp < 10 M^ to planets detectable by 
radial velocity in llda and LinI (|2005[) ^. A planet is deemed to be detectable by the radial 
velocity method if Mp > 100(a/l Af/)^/^ M® and a < 3 AU, as discussed in llda and LinI 
(12004,). The ratios are roughly 1:1.^ We scale the probability of forming a habitable planet 
by these ratios, which are binned with the following ranges of M*: 0.08-0.5, 0.5-0.8, 0.8-1.25, 
and 1.25-^1.5. The probability of forming a habitable planet with [Fe/H]>0 is given by 
Equation [B] scaled by the ratios described above. At sub-solar metallici ties, the probabilit; 



_ . ity 

of forming a habitable planet is 3% and scaled by the same method, llda and LinI (J2005I ) 



permit only one major planet around each star in their model. Therefore, a maximum of 
one habitable planet can exist around any star in our model. This study does not advance 
the notion that the stars in extra-solar systems are expected to contain a single planet. 

Given the probability of forming a habitable planet or Hot Jupiter, if a host star has 
a habitable planet and a Hot Jupiter, we assume that the habitable planet does not sur- 
vive. The probability of a star having a habitable planet without a Hot Jupiter that would 
des troy it is mod e lled a s a result of the metallicity abundance of its host star. The study 
by iHoward et al.l (|2010l ) suggests that 23% of stars harbour a short period Earth-mass 
(0.5-2 M®) planet. Although this finding does not suggest that all of these planets ex- 
ist in the HZ of their host stars, we find that ~5% of all stars host a habitable planet 
(0.1 Mq < Mp < 10 Mq) in the HZ of their host star across the disk of the Galaxy. Our 
study may underestimate the number of habitable planets that exist in the Milky Way. The 
locations of our habitable planets populated over all epochs in Models 1 and 4 are shown in 
Figure U) The inner Galaxy hosts the greatest number of habitable planets due to the high 
metallicity abundance and stellar density in the region. Note that these are the locations of 
all the habitable planets that have survived possible Hot Jupiter migrations, before deter- 
mining if they have survived SNe sterilizations. Therefore, all of these planets will not be 
suitable for advanced biological life; however, a large portion may be suitable for microbial 
life. 



3.3 Tidal Locking of Planets around Stars 

The distribution of the number and types of planets around stars of various masses is 
not fully understood at present. However, low-mass stars are not precluded from having 
Earth-mass planets. For example, the Gliese 581 system hosts super-Earths. One of these 
super-Earths exists on the outer edg e of the HZ and orbits its host star in a tidally-locked 
configuration (jvon Bloh et al.l . 120071 ). Tidal locking occurs when one side of a planet con- 
stantly faces its host star. Around low mass stars, a habitable planet must have a small 
semi- major axis to exist in the HZ. As a result of tidal friction, planets will experience tidal 
locking with their host stars on short timescales. Consequently, planets will have one side 



^Data obtained through private communication with S. Ida^ 

^The ratios are not explicitly shown as the data belongs to llda and LinI i2005l ). 
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constantly facing their host star, while the other side is in perpetual darkness. Volatiles 
on the dark-side of th e plan et may freez e out. S uch an environment is thought to make a 
planet uninhabitable ( Dolel . 1964; Kastin g et al.l . Il993). Habitability on planets that orbit 
low-mass stars has been evaluated to investigate the claim that they might be inhospitable 
to life. Research shows that heat transfer mechanisms and certain atmosph eric compo- 



sitions could preve nt volatiles from freezing out under particular conditions (|Joshi et al. 



1997[ iTarter et al.l . 120071 ). Furthermore, tidal heating could permit plate tectonics that 



woul d allow for the recycling of CO2 to prevent a runaway greenhouse effect (jBarnes et al. 
20091 ). Habitability on planets that are tidally locked with their host stars is likely to be 
different than habitability on planets that are not tidally-locked. Therefore, tidally-locked 
and non-locked planets are investigated independently in our model. 

Each star in our model is assigned an appropriate stellar mass, which is used to scale 
the probability of having a habitable planet. With regards to tidal locking, the radius 
within which an Earth-like planet would become tidall y locked is de t ermin ed, using the 



same assumptions and consequent tidal- locking line as iKasting et al.l ( 19931) . The tidal 



locking line describes the radius at which a planet will be locked into synchronous rotation 
by 4.5 Gyr for a particular stellar mass and orbital distance. Given the orbital distance of 
a planet and its host star's mass, we estima te the proba b ility o f a planet becoming tidally- 
locked using data from the simulations of llda and LinI (|2005l ). In the respective bins in 
our model, all stars with M*=0.08-0.5 host tidally-locked HZ planets; 43% of stars with 
M*=0.5-0.8 host tidally-locked HZ planets, and all stars with M*=0.8-1.5 host non-locked 
HZ planets. In other words, all M stars host tidally-locked planets, K stars host a mix of 
locked and non-locked planets and all G, and F stars host non-locked planets. 

Other factors exist that may reduce habitability around low- mass stars. Planets are 
likely to receive X-ray and UV radiation fluxes at a much higher rate around low mass s tars 



than around solar-mass stars of the same age (JLammer et al.L 120071 : ISelsis et al.l . 120071 ). In 



addition, planets that form in the habitable zones of low-mass stars may be deficient in 
volatil es as a result of high temperatures and collision rates during the planetary accretion 
phase (|Lissaueii 120071) . Furthermore, there is a tendency for low-mass stars to produce 
low-mass planets, and these planets are less likely to sustain su bstantial atmospheres and 



plate tectonic activity that is probably required for complex life (jRavmond et al.l . 120071) 



Other properties that could limit the lifespan of the biosphere are not modelled in the 
present study. For example, the HZ around a star changes as a function of time, but is not 
implemented here. Furthermore, a significant loss of radiogenic heating within the interior 
of a planet cou ld lower the CO2 levels below the threshold required for photosynthesis 
( Chaissonl . 120011 ). Hence, the biosphere may not be conducive to complex life. 



3.4 Sufficient Time for the Evolution of Land-Based Complex Life 

The only known examples of life in the universe reside on the Earth. We do not know if 
the time required for the biological evolution of complex life on Earth is representative of 
the average time required on other planets. To be able to quantify the number of habitable 
planets in the Galaxy, and compare our work to related studies, we assume that the biological 
evolution of life on other planets is identical to that of the Earth. This assumption is highly 
speculative; however, we cannot extrapolate to other evolutionary paths without knowledge 
of complex life on other planets. Our model assumes that complex life is animal life that 
dwells on land. We assume that all of these planets have land; however, if water worlds 
prove to be common, we will have overestimated the number of habitable planets. 
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The impact of a SN on life that inhabits water is not fuhy known; however, the probabihty 
of any catastrophe such as asteroid impacts, gamma ray bursts etc. entirely destroying all 
microbial life is very low. On the other hand, a SN explosion in close proximity to the 
Earth would have the capability of sterilizing all land-based complex life. This subsection 
describes how we model the timescales for the emergence of complex life, while accounting 
for the impacts of SN sterilizations. 

3.4.1 Constructing Earth- Analogs 

For complex life to exist on a habitable planet in our model, additional conditions must be 
met. These conditions are based on the milestones and consequent timescales that occurred 
on Earth and are outlined in Tabled) 

The buildup of oxygen in Earth's atmosphere coincided with th e rise of multicellular life, 
and reached present day levels between 500 and 1000 Mya ( McKavl . ll998 ). Fossil evidence of 



animal life occurs a t ^600 Mya, and mo lecular clock dates of the diversification of animals 



occurs at ~1 Gya (Catling et al.l . 120051 ) . Given that the age of the Earth is ~4.55 Gyr 
(|Allegre et al.l . Il995[ ). we assume that the rise of animal life occurred ~4 Gyr after the 
planet's formation. 

The destruction of ozone in the atmosphere of a planet caused by a nearby SN exposes 
land-base d complex life to its host star. The Earth's ozone layer was sufficiently developed 
~2.3 Gya ( Gatling et al.l . 120051) : therefore, if an Earth-analog is bombarded by the flux from 



a SN before this time, the event has no effect on the development of complex life. The 
emergence of animal life is assumed to be correlated with the rise of oxygen. Therefore, 1.55 
Gyr of continuous ozone is required for the emergence of complex life. If a SN occurs in the 
simulation between the formation of the ozone layer and the formation of animal life, we 
assume that the planet does not develop animal life until 1.55 Gyr of time elapses without 
interruption from ozone depletion events^. We require our model to reproduce major events 
that transpired on Earth to avoid speculation about life on extra-solar planets. 

As mentioned above, we ignore the production of stars in clusters. SNII that would 
normally occur within a cluster and sterilize those young planets that have formed around 
stars in the cluster do not have a significant impact on habitability. These SNII have a benign 
effect on habitability, as they occur well before ozone can form on our Earth- analogs. Those 
stars and subsequent planets nearby a SNII that are not part of the cluster that contains 
the SNII are of a distribution of ages and therefore are expected to experience sterilizing 
events. Planets orbiting these stars may have had sufficient time for the buildup of ozone 
and rise of complex life. We believe it is reasonable to assume that the background mean 
density of stars near the SNII that are not part of the cluster is a fair representation of stars 
we expect to be sterilized by the SN. This is a realistic assumption as star clusters dissipate 
on short timescales, as their stars become well mixed within the disk. We expect that an 
excess of stars above the mean stellar density at a given position as a result of clustering is 
superfluous, as this population of planets is not negatively affected by the SNII due to their 
young age. Therefore, we believe that we capture the relevant SNe sterilizations despite the 
fact that we do not model clusters. 



^For example: The ozone layer forms at 2.45 Gyr on each Earth-analog. The Earth-analog is sterilized 
by a SN at T=2.7 Gyr. The planet is not considered habitable until 4.25 Gyr in the simulation. 
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3.4.2 Timescales for the Recurrence of Land-Based Complex Life 

A planet is considered habitable if it satisfies the conditions necessary for complex life as 
described in § 13.4.11 After these conditions are met, we monitor the planet to see if they 
are maintained. If a sterilizing event occurs, we evaluate if land-based complex life can 
re-emerge. For a planet to be become habitable a successive time, the reconstruction of the 
ozone layer and the evolution of animal life must occur. 

The rise in oxygen on the Earth occurred at ^2.3 Gya, and evidence of cyanobacteria 
occurs at ^2.1 Gya. Therefore, the best-case scenario assigned for the reconstruction of the 
ozone layer is 0.4 Gyr, if cyanobacteria populations do not become extinct after a sterilizing 
event. The worst-case scenario for the reconstruction of the ozone layer is assumed to 
be 2.25 Gyr, where the ozone layer builds up along the same timescale that it took from 
the formation of Earth (4.55 Gya) to the first wave of large-scale oxygenation at 2.3 Gya. 
Therefore, the time period required for reconstituting ozone in the atmosphere is determined 
by choosing a value uniformly in the range [0.4,2.25] Gyr. 

The restoration of the ozone layer is assumed to be a prerequisite for land-based animal 
life. We assume that the best-case scenario for the redevelopment of animal life is Gyr, 
as complex life might survive in the oceans and immediately inhabit land when sufficient 
ozone is present. The worst-case scenario is 1.55 Gyr, wherein a continuous ozone layer is 
required before the redevelopment of animals. Therefore, the duration of time that elapses 
before animal life evolves and inhabits land is determined by choosing a value uniformly in 
the range [0,1.55] Gyr. 

It is not clear how long it would take for the reconstruction of the ozone layer and con- 
sequent colonization of land by complex life on exoplanets, and we can only make estimates 
of these values in the context of Earth. The timescales for the recurrence of land-based 
complex life described here are therefore more speculative than other observable properties 
that are considered in the model. 

4 Results 

We have described a model of the Milky Way which reproduces well established physical 
properties of the Galaxy and have produced a set of criteria that allows us to predict 
where and when complex life might emerge. The parameters and criteria are summarized 
in Table [3] We now present the results from our simulations. 

The results are organized as follows, ij 14. II describes the influence of SN sterilizations on 
habitability in the Galaxy without considering planet formation. § 14.21 includes constraints 
on the time required for complex life to evolve, without considering planet formation. In 
§ 14.31 we apply the SNe sterilizations and timescales for complex life to planet formation. 
§ 14.41 discusses the GHZ in terms of radial distance and height above the midplane. § 14.51 
discusses the impact stellar kinematics would have on our models. In § 14.61 we examine the 
impact of the IMF on habitability. Finally, S 14. 71 compares our results with relevant studies. 



4.1 The Effect of Supernovae Sterilizations on Habitability 

In this section, stars that are close to SN events are investigated and planet formation is 
ignored. We find that the majority of stars in our Galaxy will be bathed in flux by a 
nearby SN event during their lifetimes. Figure [5] shows the number of stars in the Galaxy 
that remain unsterilized. In this context, an unsterilized star refers to a star that has not 
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been sufficiently close to a SN event in its lifetime. Across all models, '^27%-36% of all 
stars in the Galaxy remain unsterilized. The lowest fraction of unsterilized stars in all 
models is located at R~2.5 kpc due to the high stellar density in the region. We predict 
that this fraction will be even lower at R<2.5 kpc, although we do not model this area. 
To observe how sensitive stars are to sterilizations, we double the sterilization volume of 
SNe in the simulation. When the sterilization volume is doubled, ~23% of all stars in the 
Galaxy remain unsterilized (Figure [5]- Upper panel) . Doubling the sterilization volume does 
not double the number of stars that experience a sterilization during their lifetimes because 
some stars are sterilized multiple times. This analysis indicates that a major increase in 
sterilization volume is unlikely to significantly reduce habitability within the Galaxy. 

The probability that a star remains unsterilized by a SN event ranges from '^10-90% 
in terms of galactocentric distance, and ~10-95% in terms of distance above the midplane 
(Figure [6]) across all models. We can conclude from this finding that the fraction of planets 
expected to survive a SN event is much lower towards the centre of the Galaxy, and within 
and close to the midplane. Comparing Models 1 and 4, which set the upper and lower limits 
of non-sterilized stars at R=8 kpc, we find that between ^53% and '^36% of all stars remain 
unsterilized at this radial position. 

SNII and SNIa are implemented separately to refiect the differences in their sterilization 
distances and progenitors. We find that in Model 1 SNII are responsible for more steriliza- 
tions than their SNIa counterparts. There are ~10.3x more SNII than SNIa; however, SNII 
only lead to ~1.8x more sterilizations. We estimate that an individual SNIa is ~5.7x more 
lethal than a SNII on average. In Model 4, there are ~11.7x more SNII than SNIa; however, 
SNII only lead to ~2.1x more sterilizations. An individual SNIa is roughly 5.6 x more lethal 
than a SNII in Model 4. Comparing the overall number of SNe, there are approximately 
2x more SNII and 1.8 x more SNIa in Model 4 than Model 1. The Kroupa IMF in Model 4 
produces a much higher number density of SNe, despite the lower number of stars in the 
model in comparison to Model 1. 

4.2 Timescales for Complex Life Around Main Sequence Stars 

Disregarding planet formation, we add a 4 Gyr time constraint corresponding to the timescale 
required for the emergence of complex life. In this subsection we analyze the positions of 
stars that are suitable for complex life when considering SNe sterilizations. Hence, the star 
must have 1) a main sequence lifetime of at least 4 Gyr, and 2) remain unsterilized for a 4 
Gyr period. If a star is sterilized during this period, the biospheres of orbiting planets could 
be affected. In this scenario, SNe will have the effect of delaying the emergence of complex 
life. The 4 Gyr constraint on those stars suitable for habitable planets is relaxed when 
planet formation is considered in Section 14.31 where the assumptions for the development 
of complex life in Section [33] are adopted. 

Figure [7] presents the fraction of stars that are unsterilized for at least one period of 4 
Gyr. Interestingly, the inner Galaxy contains a fraction of sterilized stars that nevertheless 
exhibit a 4 Gyr unsterilized period that is suitable for complex life. Across all models, the 
lowest fraction of stars that are 4 Gyr of age and remain unsterilized are located in the 
inner Galaxy (Figure [T]- Upper panel). This is expected, as the majority of stars in the inner 
Galaxy are sterilized (Figures [5] and [6]) . The reverse phenomenon is demonstrated in the 
outskirts of the Galaxy in all models; the majority of stars suitable for complex life are 
unsterilized and at least 4 Gyr of age. Due to the declining stellar density as a function of 
radial distance, there are more stars in total that would permit the emergence of complex 
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life in the inner Galaxy than the outer Galaxy. A high stellar density corresponding to a 
high SN rate would completely inhibit the emergence of complex life. At R>2.5 kpc there 
are no locations where this occurs, given the assumptions made in our model. 

In the middle and bottom panels of Figure [71 the fraction of stars with 4 Gyr lifespans 
declines at R^12 kpc, as many of the stars in the region are not 4 Gyr of age; therefore, 
these stars are less habitable than those at R<12kpc. Furthermore, the fraction of stars that 
have 4 Gyr periods where they are unsterilized by SNe across the disk is higher in Model 1 
in comparison to Model 4, as the Salpeter IMF produces fewer SNe than the Kroupa IMF. 

4.3 The Effect of Metallicity on Planet Formation With Respect 
to Habitabihty 

We now present our model of habitabihty with regards to planet formation as influenced by 
the metallicity gradient in the disk of the Milky Way galaxy. Section H?^ discussed stars that 
are good candidates for habitable planets. Two criteria in assessing the habitabihty of these 
stars were implemented: 1) the star must have a main sequence lifetime of at least 4 Gyr 
and 2) the star must remain free from sterilizing SNe events for a period of 4 Gyr. Requiring 
criteria (2) is a strict prerequisite and does not reflect the impact SNe would have on our 
Earth-analogs given the major events that have transpired in Earth's history. For example, 
a SN would have no effect on a planet that has not formed an ozone layer. These milestones 
were reviewed in § 13.41 and the results of their implementation are presented in this section. 
In this section, we highlight the results of Models 1 and 4 in parts of our analysis, as they 
most closely constrain the local number density in the solar neighbourhood. Note that the 
Galaxy is best represented by the Kroupa IMF, and therefore we can expect Model 4 to 
best represent habitabihty within the Galaxy. 

We find that the greatest number of habitable planets exist in the inner Galaxy in all 
of our models. More specifically, 50% of the habitable planets lie at R<4.1 kpc and R<4.4 
kpc in Models 1 and 4 respectively. The high stellar number density towards the centre 
of the Galaxy and the high degree of metallicity are responsible for permitting more stars 
to have habitable planets in this region in comparison to the middle and outer Galaxy 
(Figure [8|). Moreover, star formation occurs earlier in the inner Galaxy, permitting longer 
periods for the emergence of complex life. Furthermore, when comparing the age distribution 
of stars in the inner, middle, and outer regions of the Galaxy, the inner region allows for 
a greater chance of attaining habitable conditions in comparison to younger stars in the 
middle and outer regions of the disk. The high stellar number density in the inner Galaxy 
and subsequent sterilization rate does not completely frustrate the emergence of complex 
life in t he region. As mentio ned previously, only one major planet forms around a star in the 



work of llda and LinI (J2005l) : therefore, the number of habitable planets is underestimated 
here. 

With respect to the fraction of habitable planets, we predict that between s»1.2% 
(Model 4) and wl.7% (Model 1) of all stars host a habitable planet. More specifically, in 
Model 1, 1.60% of stars host a tidally-locked HZ planet and 0.09% of stars host a non-locked 
HZ planet. In Model 4, 0.9% of planets are tidally locked, and 0.3% are in a non-locked 
configuration with their host star. The Salpeter IMF (Model 1) leads to more low mass stars 
whose companions become tidally locked on short timescales. However, the Kroupa IMF 
(Model 4) leads to higher mass stars, permitting more planets in a non-locked configuration. 
In our model, the mix of locked and non-locked planets is constant throughout the Galactic 
disk. The area with the greatest fraction of habitable planets over all epochs is located 
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towards the centre of the Galaxy, as shown in Figure [SI 

From Sections 14.11 and 14. 2[ we see that SN sterihzations on their own make the inner 
Galaxy the least hospitable for complex life. However, regarding the planet-metallicity 
correlation without the effects of SNe makes the inner Galaxy the most hospitable for 
complex life. When both factors are taken into account, the inner Galaxy is ~10x more 
hospitable than the outer Galaxy (Figure IH]). This finding indicates that the impact of 
metallicity on planet formation appears to dominate over the effects of SN sterilizations. 
Furthermore, the inside-out scenario of Galaxy formation permits the inner region to be 
more habitable than the outskirts. Neither SN sterilizations nor metal-poor environments 
are capable of rendering any region inhospitable to complex life at the present day. 

With respect to time, we investigate 1) the formation date (birth date) of habitable 
planets, and 2) the time in a planet's history when it becomes habitable. Figure [TOl shows 
the locations and birth dates of planets that are habitable at the present day. All of our 
models indicate that most of the habitable planets are located in the inner Galaxy. We find 
that our solar system is younger and distant from the densest regions that host habitable 
planets. In addition, our results indicate that the average birth date of planets occurs 4.06 
Gya, and 3.94 Gya in Model 1 and 4 respectively. Thus, many planets are too young to 
permit the evolution of complex life. Over the next Gyr, many more planets are expected 
to attain habitable conditions given the assumptions made in our model. 

With respect to when planets attain habitable conditions, the top panel of Figures [TT] 
and [13 demonstrate that habitable planets have emerged for the past ~6 Gyr near the 
Galactic centre (R>2.5 kpc), whereas they have become habitable at ^A Gya in the middle 
region, and at ~l-2 Gya in the periphery. The SFR experienced in the last few billion 
years, coupled with increasing levels of metallicity, suggests that many more planets will be 
conducive to complex life in the future. 

Tracing the habitability history of planets over all epochs shows that at the present 
day, the fraction of stars hosting planets with habitable conditions is greatest at R~2.5 kpc 
(Figures [TT] and [T21- Top Right panel). The inner Galaxy has the greatest fraction of stars 
with a habitable planet integrated over all periods of time (Figure [9]) , and at the present 
day. 

4.4 The Location of the GHZ 

The GHZ is thought to be affected by an inner boundary that is determined by hazards 
to planetary biospheres, and a n outer boundary set b y the minimum amount of metallicity 
required for planet formation ( Gonzalez et al.l . 120011) . This outlook describes the GHZ as 



a function of radial distance. However, we model concentrations of habitable planets with 
respect to galactocentric distance (R) and height above the midplane (z). 

The bottom left panels of Figures [TT] and [12] illustrate that the greatest number of 
habitable planets over all epochs is located in the inner Galaxy, within and surrounding the 
midplane. This region coincides with the greatest stellar number density in the model. Most 
planets in this densely populated area will be sufficiently close to SNe, but the frequency 
of ozone depletion events is not high enough to permit large volumes of the Galaxy to be 
sterile for long timescales. 

In the bottom right panels of Figures [TT] and [12] it is observed that the region with the 
greatest fraction of stars with habitable planets over all epochs is located in the inner Galaxy, 
well above the midplane, centred at a height of z~1.5 kpc. The same radial mix of metals 
at this location above the midplane combined with the low stellar density in comparison 
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to the density found closer to the midplane suggests that SNe have a significant impact on 
habitabihty. The fraction of stars with habitable planets above the midplane is a factor of 
a few greater than the fraction at the midplane at R~2.5 kpc. 

We observe that the greatest number density and fraction of habitable planets is located 
in the inner Galaxy, integrated over all epochs. These findings suggest that our location 
in the Galaxy is not particularly favourable under this GHZ paradigm. We find that the 
GHZ, defined as the position with the greatest number of complex life supporting habitable 
planets, is located in the inner Galaxy, within and surrounding the midplane. 

4.5 Stellar Kinematics 

The motions of stars are not considered in our model. Solar systems passing through dense 
regions of the Galaxy are more likely to be sterilized by SNe. The lack of stellar kinematics 
would be a major concern if we found a continuously sterilized zone (CSZ), wherein an entire 
section of the Galaxy is uninhabitable as a result of transient radiation events, effectively 
sterilizing stars that pass through the region. However, the fraction of habitable planets in 
the densest region (at the midplane in the inner Galaxy) is not significantly lower than most 
of the other parts of the Galaxy. Therefore, if the orbits of stars were modelled in the present 
study, we would not expect a major decrease in habitable systems as a result of traversal 
through high density regions. Given the results presented in this study, neither radial or 
vertical motion presents a major concern. In the radial case, if a star in the midplane (z^O 
pc) moves inwards towards the Galactic centre, the star will be in an area where the fraction 
of stars with habitable planets increases (Figures [TT] and [T^Bottom left panels), therefore, 
the SNe rate cannot be great enough to inhibit habitabihty in the majority of solar systems 
with radial components to their motion. In the vertical case, if a star in a low density region 
above the midplane (e.g. R^2.5 kpc and z~2 kpc) plunges towards the midplane, the 
fraction of habitable planets only decreases by a factor of '^2-3 (Figures [Til and \T2\ Bottom 
right panels). The slightly decreasing fraction of habitable planets towards the midplane 
at R~2.5 kpc suggests that vertical motion cannot be responsible for a major decrease in 
habitabihty. If vertical stellar motions are taken into account, the vertical gradient of the 
fraction of habitable planets is likely to be reduced. Further research on habitabihty in 
the Galactic bulge in combination with stellar kinematics is necessary to determine if dense 
regions have a significant impact on habitabihty. If a large CSZ exists in the inner disk or 
the bulge, stars that reside in the region, or those stars in eccentric orbits that pass through 
the region, may not have planets that can support complex life as defined in this study. 

4.6 The Impact of the IMF on Habitabihty 

We present four models to see how sensitive our model is to varying properties of the Galaxy. 
A Salpeter IMF is utilized to assign masses to the stars in Models 1 and 3 and the Kroupa 
IMF is used in Models 2 and 4. 

Comparing Models 1 and 4, we see that the former has a lower SNe rate than the latter. 
There are ~1.8x more SNII and '^1.7x more SNIa in the models with the Kroupa IMF 
in comparison to those with the Salpeter IMF. Taking the SNe rates in consideration on 
their own makes the Galaxy less hospitable for complex life when the Kroupa IMF is used 
to assign masses to our stars. However, when comparing the fraction of stars that host 
habitable planets in a non- locked configuration (Figure [9]) , the Kroupa IMF produces more 
high mass stars permitting many more planets to orbit their hosts in a non-tidally locked 
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configuration in their HZ in comparison to planets formed in the models containing the 
Salpeter IMF. Evidence suggests that the Galaxy follows the Kroupa IMF; therefore, we 
expect '^0.3% of stars to host a planet capable of harbouring complex life that orbits in a 
non-locked configuration. 



4.7 Comparison with other Studies 



PrantzosI (|2008l ) suggests that the concept of the GHZ may have very little significance, 
and that it should be considered only as a framework to organize our ideas about life 
in our Galaxy. While it is understood that further research is required to better answer 
questions related to habitability on the Galactic scale, we do believe that these factors 
can be quantified at present with useful results. It is certain that a sufficiently nearby SN 
would have a sterilizing impact on land-based life on Earth, and the evolution of biologically 
complex life took a substantial amount of time on our planet. Studies of this nature imply 
uncertainties that require us to make reasoned assumptions concerning habitability. In the 
future, these assumptions will be replaced with observational fact. 

In this subsection we compare key differences between the present study and the related 
literature. Unlike previous approaches, our simulation models each star individually. While 
such a star-based model incurs increased memory and computational costs, it potentially 
affords more realistic and informative results than traditional probabilistic models that 
aggregate over stellar populations. The contrasting methods and results are reviewed, some 
of which have been described or alluded to in previous sections. 



4.7.1 Model Milky Way Galaxies 



The buildup of the Galactic stellar mass rep orted in iPrantzos (120081) and the present study 
is roughly consistent. The SFH employed in IPrantzosI (|2008[ ) . iLineweaver et al.l (|2004l ). and 
our model, has an early burst of star formation as illustrated in the upper and middle 
panels of Figure [TJ The average metallicity assigned to each sta r in our model (F igure [1] 
- lower panel) , is roughly consi stent with the rn etallicity found in ILineweaver et al.l (J2004l ) . 
The measure of metallicity in IPrantzosI ( 20081 ) is different than ours; therefore, a direct 
comparison is not possible. 



4.7.2 Supernovae Rates and Sterilization Distances 

Our research addresses the dang er effect of SNe using a different approach than that of 
IPrantzos (2008) and Lineweaver et al. (2004) . Numer ous studies in the field show that SNIa 



and SNII have different distributions of luminosities ( Richardson et al.1 . 120021 : iWang et al. 



20061 ). Our model accounts f or this usin g a distribution of sterilization distances for SNII and 



SNIa. Conversely, |Prantzos| ([2008 ) and ILineweaver et al.l ( 2004) assess the SN danger factor 
as a time integrated supernova rate. Treating SN events in our model on an individual basis 
within the pre-existing stellar population improves on the statistical approaches of other 
studies, as it is a more realistic, less Earth-centric method. Our model does not assess the 
danger posed by SNe as normalized to the Earth's radial position. 



4.7.3 Metallicity and Planet Formation 

There are few studies that predict t he number of habitabl e planets in the Milky Way 
(JLineweaven . 1200 ll iBounamam et al.l . 120071 : IGuo et al.l . 120091 ). We use a metallicity pro- 
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file (JNaab and Ostrikeii 120061 ) to assign a metallicity to ea ch star in the m odel. We use this 
metallicity to model planet formation. Considering that IPrantzoa (|200 8^ uses a constant 
probability of stars forming an Earth-like planet (40%), as influenced by [Lincwcavcr (2001), 
our model contains far fewer habitable planets. Moreo ver, the danger of a Hot Jupiter 
inhibiting an Earth-like planet in our study is similar to IPrantzod ( 20081 ). as we share the 
results by F ischer and Valcnti (2005 ), producing a much lower Hot Jupiter danger effect 
than that of iLineweaver et al.l (|2004l) . 



4.7.4 Comparison of Results 

The GHZ is a n area in the Galaxy that contains stars with the highest potential to harbour 
complex life. ILineweaver et al.l (J2004l) find the GHZ to be an annular region between 7 
and 9 kpc at the present day. Our results fundamentally disagree with ILineweaver et al.l 



(|2004[ ). as we find that the greatest number of habitable pl anets exist in the in n er Ga laxy. 
However, we recognize that the mod el Galaxy produced bv lFenner and GibsonI (|2003[ ) and 
employed in Line w eaver et al.l (|2004[ ) is more advanced than our model in certain regards. 
Their model contains a more accurate depiction of the observed properties of the Galaxy 
such as considering the SFR in conjunction with spiral arm motions. 

We expect that as the metallicity increases, the entire disk is expected t o harbor a greater 
number and fraction of habitable planets. This result is similar to that of lPrantzoa ([2008) ■ 
Moreover, we find that the number of habi table planets varies as a fu nction of height above 
the midplane. This was not investigated in ILineweaver et al.l (2004j or Prantzos| ( 20081) . and 
th erefore we observe that t he morphology of the GHZ is not a n annular r egion as suggested 



Lineweaver et al.l (|2004l ). or the entire disk as discussed in IPrantzoJ (|2008l ). but rather 



consists of a region highly dependant on radial distance, located near the centre of the 
Galaxy. 

Habitability may be reduced at R^2.5 kpc as a result of the high stellar density in 
combination with stellar motions, and will be investigated in future research concerning 
the Galactic bulge. Nonetheless, our prediction that many habitable planets exist at the 
midplane in the inner Galaxy suggests that habitable conditions are possible in high density 
regions. 

The dominant paradigm of habitability on the Galactic scale is certain to change in this 
field as we learn more about our Galaxy and habitable planets. Studies relating to the 
distribution of habitable planets, extraterrestrial life and colonization of the Gala xy (see 
iBouuamam et al. , 2007; Forgan, 2009) currently rely on the canonical model bv Lineweaver et al 
(J2004l ). Our revised conception of the GHZ has implications for these and other related stud- 



5 Conclusions 

We present a model of habitability within the Milky Way to predict the region(s) in the 
Galaxy that are expected to favour the emergence of complex life. We do not find that 
the inner Galaxy is entirely inhospitable to life; in fact, the greatest number of habitable 
planets are found in this region. The metallicity in the inner Galaxy produces a high 
planet formation rate for long timescales that dominates the negative impact SNe have on 
habitability. Furthermore, we observe that over all epochs in our favoured Model (Model 4) , 
1.2% of all stars in the Galaxy host a habitable planet (including both tidally locked and non- 
locked configurations). The fraction of stars with a habitable planet ranges from ~0.25% 
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in the outer Galaxy to ^2.1% in the inner Galaxy. Considering that we find the greatest 
number of habitable planets to be located in and around the midplane at R~2.5 kpc, and 
that the greatest fraction of stars that support habitable planets exist above and below the 
midplane at this radial position, we find that the inner region of the Galaxy may support 
the greatest number of planets conducive to complex life. More specifically, our findings 
are at odds with the notion that the GHZ is s haped like an annular ring at R~7-9 kpc at 
the present day, as found by iLineweaver et al.l (|2004l ). Similarly to IPrantzosI (|2008l) . we find 



the greatest number of habitable planets to exist in the inner Galaxy, with the exception 
that we observe habitable planets at this radial position to be strongly dependent on height 
above and below the Galactic midplane. 

Given that we find the greatest concentration of habitable planets to be located in the 
inner Galaxy, further research on the disk and overlapping bulge at R<2.5 kpc is war- 
ranted. The inner bo undary of the GHZ is defined by hazards to a planet's biosphere 
( Gonzalez et al.l . 120011 ): however, we find that an inner boundary does not exist when mod- 



elling the Galactic disk at R^2.5 kpc. By modelling the region at R^2.5 kpc, we might find 
an inner boundary, or we might not find a boundary at all. 

Studies of habitability on the Galactic scale will imp rove in the future as Earth- like plan- 
ets are detected from studies such as the Kepler mission ( Borucki et al.l . |2008|) , and will bene- 



fit fro m an improved understanding of our Galaxy from the GAIA mission (|Perrvman et al. 



20011 ). and others. 

The results of planet finding missions will yield estimates of the total number of planets 
in the HZ of their host stars across the Galactic disk. The total number of habitable planets 
estimated here is not the focus of our research; rather, we highlight and predict that these 
planets exist preferentially in the inner Galaxy, and that they have the ability to survive 
SN sterilization events for periods conducive to the rise of biologically complex life. 
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Figure 1: Temporal profile of the SFR and metallicity in our model at six galactocentric 
regions: 2.5, 5, 7.5, 10, 12.5, 15 kpc ordered from top to bottom. Upper panel: Temporal 
profile of the SFR in Model 1. The curves indicate a burst of star formation in the inner 
Galaxy that declines with time. The SFR at 7.5 kpc had a burst of star formation early 
in the history of the Galaxy and flattens within the last few Gyr. Middle panel: Temporal 
profile of the SFR in Model 4. The SFR is similar to that of Model 1, except the SFR 
is slightly lower in the inner Galaxy, and higher in the outskirts. Lower panel: Temporal 
profile of the average metallicity in our models. The inner Galaxy produces metals early as 
a result of the burst of star formation in the region, whereas the build up of metals takes 
longer in the outer Galaxy. 
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Figure 2: The present day distribution of mctallicity for all stars in our model that have 
formed at the sola. r radiu s (R'^8 kpc). This distribution is consistent with the model of 
Naab and Ostrikeij (|2006l ). 




Figure 3: The distribution of sterilization distances for SNII (circle) and SNIa (square). The 
curves represent the distribution of sterilization distances, corresponding to the distribution 
of absolute magnitudes of SNII and SNIa. The x value is used to select a sterilization dis- 
tance for a given SNII or SNIa, where a; is a random number generated between and 1 
in the Monte-Carlo simulation. It is not surprising that the range of sterilization distances 
for SNIa is less than that of SNII, reflecting the view that SNIa are considered to be stan- 
dard candles. SNIa are more luminous on average and therefore have a greater average 
sterilization distance than SNII. 
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Figure 4: The positions of all of the habitable planets formed over all epochs in Model 1 
(upper curve) and Model 4 (lower curve). These planets have 0.1 M^ < Mp < 10 Mq, exist 
in the HZ of their host stars and have survived the negative effects posed by Hot Jupiters. 
The discrepancy between the number of planets formed in Models 1 and 4 is a function of 
the IMF, where the Salpeter IMF in Model 1 produces more stars, and hence, more planets 
than the Kroupa IMF utilized in Model 4. This may give us a rough idea of the distribution 
of planets capable of hosting microbial life. 
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Figure 5: The number of stars in the Galaxy that would not have their respective solar 
systems sterilized by a SN is plotted as a function of galactocentric distance for each model. 
Upper panel: Models 1, la and lb are plotted. The top curve represents the total stars in 
the model. The curve with the square markers represents those non-sterilized solar systems 
with our default sterilization distance. We find that ~33.3% of stars are never sterilized 
by a SN event. Stars in the inner Galaxy are more likely to be sterilized than those in the 
outskirts. Our sensitivity analysis indicates that when the sterilization volume is halved, 
^45% of all stars in the Galaxy remain unsterilized (circle markers-Model la). When we 
double the sterilization volume (triangle markers- Model lb), ^^23% of all stars in the Galaxy 
remain unsterilized. Middle panel: Model 3 is plotted. The top curve represents the total 
stars in the model. In this model ~36% of all stars remain unsterilized (bottom curve). 
Bottom panel: The total stars and unsterilized stars for Model 2 (square markers) and 
Model 4 (circle markers) are plotted. In Model 2, ~27% of all stars remain unsterilized, 
whereas in Model 4, ^29% of all stars remain unsterilized. 
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Figure 6: Upper panel: The fraction of stars that are not sterilized by a SN event as a 
function of galactocentric distance. Models 1-4 are labeled with the square, circle, triangle, 
and diamond markers respectively. Only a small fraction of stars remain unsterilized by a 
SN event in the inner Galaxy, whereas almost all stars are unaffected in the outer Galaxy. 
All models converge in the inner Galaxy where ~10% of all stars are not sterilized. Likewise, 
in the outer Galaxy, all models converge, yielding >90% of all stars to remain unsterilized. 
At R:=8 kpc, between ~36% (Model 4) and ^^53% (Model 1) of all stars remain unsterilized 
at this radius. Lower panel: The fraction of stars that are not sterilized by a SN event as 
a function of height above the Galactic midplane. Only a small fraction of stars remain 
unsterilized by a SN event at the midplane, whereas almost all stars remain unsterilized at 
z~2.5 kpc. 
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Figure 7: Upper panel: The total fraction of stars that have 4 Gyr periods where they 
are unsterilized by SNe is plotted versus radial distance. The square, circle, triangle, and 
diamond markers correspond to Models 1, 2, 3, and 4 respectively. The inner Galaxy has the 
lowest fraction of stars suitable for complex life. At R=8 kpc, between ~15% (Model 4) and 
^^30% (Model 1) of stars have 4 Gyr periods where they are not sterilized by SNe. Middle 
panel: Model 1 is plotted. The curve with triangle markers represents the fraction of stars 
that are never sterilized; the curve with circle markers represents the fraction of stars that 
have been sterilized, but still have a period of 4 Gyr absent from sterilizing events; and the 
curve with square markers (triangle + circle) represents the total fraction of habitable stars, 
which have 4 Gyr periods unsterilized by SNe and is also shown in the upper panel. The 
majority of stars suitable for complex life in the inner Galaxy have been sterilized. Bottom 
panel: The same as the middle panel; however. Model 4 is plotted. We find that roughly 
half of the stars that are capable of supporting complex life are never sterilized at 8 kpc. 
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Figure 8: The number of habitable planets is plotted versus radial distance. Upper panel: 
Models 1 (square markers) and 3 (circle markers) are plotted (models with a Salpeter IMF). 
The curves indicate the total number of habitable planets that occur over all epochs of 
time, including both tidally-locked and non-locked configurations. The greatest number 
of habitable planets is located in the inner Galaxy. In Model 1 «1.7% of all stars host 
a habitable planet, 1.60% of which are considered to be tidally locked to their respective 
host star and 0.09% are in a non-tidally locked configuration. We find that 50% of the 
habitable planets in Model 1 lie at R<4.1 kpc. Model 3 is similar to that of Model 1. Lower 
panel: Models 2 (square markers) and 4 (circle markers) are plotted (models with a Kroupa 
IMF). The results are very similar to those of Models 1 and 3, where the greatest number of 
habitable planets is located in the inner Galaxy. In Model 4 we find that «1.2% of all stars 
host a habitable planet (0.9% are tidally locked, and 0.3% are in a non-locked configuration 
with their host star). 
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Figure 9: The fraction of stars with a habitable planet over all epochs in the Galaxy is 
plotted as a function of radial distance. Model 1 is plotted in the upper panel, and Model 4 
in the lower panel. The curves from top to bottom in each panel indicate the total number 
of habitable planets, those in a tidally-locked configuration, and those not tidally-locked to 
their host star. There is a higher probability of a star having a habitable planet in the inner 
region of the disk. The probability declines in the outer Galaxy, as the metallicity is not 
sufficient to produce high planet formation rates, despite the much lower SN sterilization 
rate. Furthermore, planets form at an earlier time in the inner Galaxy than in the outskirts, 
allowing for more planets to attain habitable conditions, despite the higher SNe rate. The 
distribution of stellar masses and subsequent probabilities of planet formation around stars 
of different masses indicates that the majority of habitable planets will be in a tidally-locked 
configuration. From the plots, the IMF has a major influence on the fraction of planets that 
are not tidally locked to their host star. 



33 




Radial Distance (i^pc) 
Model 3 





10 12 14 



Radiai Distance (kpc) 
IVlodel 4 



6 8 10 12 

Radiai Distance {i<pc) 




6 8 10 12 

Radiai Distance (kpc) 



Figure 10: The number of habitable planets per pc is plotted as a function of radial distance 
and birth date for each model. We trace the history of each habitable planet to determine 
in which periods they remain habitable, and plot those that are habitable at the present 
day. It is clear, given the assumptions made in our models, that habitable planets are 
most prevalent in the inner Galaxy. Therefore, the inner Galaxy is expected to favour the 
emergence of complex life and the region is expanding outwards with time. The majority of 
stars that host habitable planets in our model arc expected to be older than the Sun. This 
figure can be directly compared with Fig. 3 in ILineweaver et al.l ( 20041 ). 
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Figure 11: In this Figure, Model 1 is plotted. Top Left panel: The number of planets 
that are habitable (tidally-locked and non-locked) as a function of radial distance and time. 
We trace the history of each habitable planet to determine in which periods they remain 
habitable. It is clear, given the assumptions made in our model, that the number of habitable 
planets is greatest in the inner Galaxy, at all epochs. Top right panel: The fraction of 
stars with a habitable planet as a function of time and radial distance. We trace the 
history of each habitable planet to determine in which periods they remain habitable. Given 
the assumptions made in the model, the inner region of the Galaxy exhibits the highest 
probability of having habitable planets at the present day. At R^5-ll kpc, at the present 
day, the entire range has roughly the same probability of having habitable planets. Lower 
Left panel: The number of habitable planets integrated over all epochs as a function of 
radial distance and height above the midplane. We predict that the position in the Galaxy 
with the greatest number of habitable planets is located in and around the midplane in 
the inner Galaxy. Lower Right panel: The fraction of habitable planets integrated over all 
epochs as a function of radial distance and height above the midplane. The region with the 
greatest fraction of habitable planets exists above the midplane in the inner Galaxy. The 
high metallicity that produces a high planet formation rate and the lower stellar density that 
exists well above the midplane at this radial position permits a greater fraction of habitable 
planets to form. 
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Figure 12: The same as Figure [TTl except Model 4 is shown. The trends are very similar in 
both models. 
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Table 1 


.: Descript: 


ion of the models 








Model^ 


Stellar Number Density Distribution 


IMF 


Sterilization Volume 


Number of Stars^ 


n(30 


pc, 8 kpc)^ 


1 


Carroll and Ostlie (2006) 


Salpeter 


1 


1.03 




0.155 


2 


Carroll and Ostlie (2006) 


Kroupa 


1 


0.366 




0.055 


3 


Juric et al. (2008) 


Salpeter 


1 


1.13 




0.220 


4 


Juric et al. (2008) 


Kroupa 


1 


0.402 




0.078 


la 


Carroll and Ostlie (2006) 


Salpeter 


0.5 


1.03 




0.155 


lb 


Carroll and Ostlie (2006) 


Salpeter 


2 


1.03 




0.155 



^ All models employ the metallicity profile and SFH of .Naab and Ostriken l|2006l ). The total disk mass in all of the models 

is 4.2 X 10^°Mq between kpc <R<15 kpc. 
■^ Units of lO'^^ stars between 2.5 kpc <R<15 kpc. 
^ The stellar number density at the solar neighbourhood in units of stars pc"''. 
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Table 2: Major events in Earth's History that inform the timeseales for the emergence or 
re-emergence of complex life 



Event Name 


Time of 


occurrence (Gya) 


Rise of metazoan life 




0.75 


Formation of the ozone layer 




2.3 


Evidence of cyanobacteria 




2.7 


Formation of the Earth 




4.55 
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Table 3: The main characteristics that are apphed to each star in the model 



Stellar Property 



Formula or Input Parameter 



Initial Mass Function 



Salpeter or Kroupa 

The SFH given bv lNaab and Ostrikei) (120061 )" 



Birth date 



After the main sequence lifetime expires (Equation [TJ 



Death date 



Position 



Stel lar number d ensity distribution of iCarroU and Ostlia 1^2006, ) or 

Juric et alJ (|2008ll ,_^^ 

The metallicity distribution from Naab and Ostriker' (2006) 



MetaUicity 



SNII 



All stars with main sequence M > 8 Mq (Kennicutt , 1984 ) 
1% of white dwarfs (PritchelLeLal ":2QQi) 



SNIa 



Sterilization distance 



Normali zed to t he a verage 
pc ( Gehr els et "all . l2003l 'l 



SNII sterilization distance of 8 



Having a habitable planet A function of the planet-metallicitv correlation ('Fisc her and Valentil 

|2005 ) and a model of solar system formation (Ida and Lid . |2005| ) . A 
habitable planet has 0.1 A/® < Mp < 10 M®, and exists in the HZ 
of its host star 

Tidally-locked planet Those with an appropriate mass, corresponding to the tidal-locking 

line of Kasting et al. (1993) 

Ozone reconstruction time A value chosen uniformly in the range [0.4,2.25] Gyr 

Animal re-evolution time A value chosen uniformly in the range [0,1.55] Gyr 
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